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ABSTRACT 






Acoustic waves of finite extent were propagated into a dc glow 
discharge. It was found that compression phases of these waves reduced 
light emitted by the discharge, while rarefaction phases increased this 
light. In the compression phase of the acoustic wave these changes were 
attributed to increased losses of electron energy through an increase in 
the collision frequency. Electron energy was increased in the rarefaction 
phase of the acoustic wave as collision frequency was reduced. 

Fast and slow waves of stratification were excited by the acoustic 
waves. These waves of stratification had opposite signs of light in- 
tensity variation when excited by compression and rarefaction phases of 
the acoustic wave. 
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CHAPTER I 



INTRODUCTION 

Electric and magnetic fields have been widely used for investigation 
of the properties of glow discharges because of their strong interaction 
with charged particles. In recent years, the interaction of acoustic 
waves with the discharge plasma has received some attention. These ex- 
periments have been limited in scope and little theoretical work has 
accompanied them. This work is an experimental study of the effects of 
an acoustic wave group of finite extent upon a neon glow discharge 40 cm 
long, in a 2.5 cm diameter tube in the pressure range 0.7 to 4.0 torr, 
and operated at currents up to 35 ma. 

In the discharge of electricity through low pressure noble gases, 
spontaneous oscillations of the discharge current are nearly always ob- 
served. Frequencies of these oscillations range from 0.5 to 50 kHz. 
Photomultiplier signals indicate that the light intensity of the positive 
column also oscillates with this frequency. These oscillations are 
caused by the anode to cathode motion of light areas in the discharge 
which are known as moving striations. 

Studies of moving striations have made great progress in recent 
years by investigating transients introduced into a homogeneous discharge 
through impulsive modulation of the current, or through application of a 
high voltage pulse to a ring encircling part of the positive column. 

These transients have the form of a wave group or wave packet which prop- 
agates toward the anode with a group velocity u, while the component 
waves have phase velocity v directed toward the cathode. It is therefore 
a backward wave. Individual waves have the characteristics of moving 
striations, and the group has been termed a ’’wave of stratification”. 
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Figure 1-1 shows schematically the traces observed by a photomulti- 
plier at two points in the column (1-1 (a), (b)). From these the propaga- 
tion in the direction of the anode can be seen. Part (c) of this figure 
shows development of the wave of stratification in time and space, with 
positive X chosen in the direction of the electric field. In this dia- 
gram, the dark areas correspond to the moving striations which are present 
within a wave packet. Striation motion toward the cathode is evident, 
and the striation amplitude is depicted by the width of the line. 

Under conditions in which this wave receives large amplification 
during propagation along the column, the oscillations of current which 
occur as it arrives at the anode cause excitation of additional waves of 
stratification. Further increase of the amplification causes the striations 
to be self-excited, with striations of largest amplitude appearing at the 
anode end of the discharge. 

Pekarek first reported these waves in 1954, and his ensuing theoret- 
ical and experimental works have established several important character- 
istics: 

1) the waves are ionization phenomena, dependent upon the 
change of ionization rate with electron temperature; 

2) a single slow and two fast waves exist, with the slow wave 
dependent upon stepwise ionization, the fast waves on direct ionization; 

3) the slower of the fast waves is associated with motion of 
an atomic ion, the faster with the motion of a molecular ion; 

4) any localized perturbation of ionization rate of sufficient 
amplitude will cause the wave to originate; perturbations of opposite 
sign produce waves with opposite courses of light intensity; 

5) the formation of space charges due to higher mobility of 
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(a) 



(b) 





Figure 1-1. (a) Time modulation of light intensity in anode 

region by wave of stratification. 

(b) Time modulation of light intensity nearer 

cathode. 

(c) Time-space development of light intensity 
by wave of stratification. 
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the electrons, and the consequent change of electron temperature, cause 
the wave to propagate; 

6) striations with small separation (or wavelength, ^ ) are 
damped by longitudinal ambipolar diffusion, those with long wavelength 
by relaxation of electron temperature. 

Pekarek has shown that any spatial perturbation of the steady-state 
ion density causes a localized space-charge region to be formed. j 

On this basis the formation and motion of striations can be explained by 
the following chain of events: 

a) ion density is increased; 

b) a space charge region is formed 5 

c) electric field is increased on the anode side of this region, 
decreased on the cathode side; 

d) electron temperature changes with the electric field, causing 

e) increase of light toward cathode, decrease toward anode; 

f) ions move toward the cathode under the influence of the main 
electric field. 

Propagation of the wave toward the anode can be explained by considering 
the decrease of electric field on the anode side of this disturbance. 

In this case: 

g) electron temperature and ionization rate are decreased, 

causing 

h) a reduction of light intensity and a negatively charged 



reg ion; 

i) this region of negative charge enhances the electric field 
on its anode side, causing 

j) increase of ionization rate and ion density. 
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Step (j) was the initial condition of the disturbance. Since the 
first ion density perturbation caused a second similar region to form 
on the anode side, repetition of this cycle describes propagation of the 
wave in the direction of the anode. 

Properties of the wave can be related to microphysical processes in 
the discharge through the relaxation time which is defined as 



To 




1 

U + V 



( 1 ) 



This relaxation time is characteristic of the slowest process in the 
chain of events which form and propagate the wave. 

In the case of the fast wave, it is the time required for the excess 
space charge to vanish through diffusion losses, while in the slow wave 
it is the sum of this diffusion time and the time required to ionize the 
excess of ions produced by the increase of electron temperature. 

Dispersion of this wave in the inert gases is such that the higher 
frequency components of the packet generally appear earliest at a given 
position. In time-space photographs, this type of dispersion causes the 
striation paths to be curved convex with respect to origin. 

Amplification of the wave is dependent upon the change of ionization 
rate with temperature, and upon the spatial phase difference between the 
course of electric field variation and that of the electron temperature. 

[-] 

Thus the essential properties of the wave of stratification are 
well understood. Detailed mathematical formulation of the wave of strati- 
fication has been carried out by Pekarek and various co-workers ^41 - 60 j, 
and the theory for the case of zero Debye length has been found to des- 
cribe this case reasonably well. However, all perturbation work has 
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been accomplished electrically with the exception of some work by Sicha 
et. al., who used square-wave modulation of the power input to re-entrant 
RF cavities to excite the wave of stratification, thus directly per- 
turbing the ion density locally, Other areas of the column are 

affected, however, as standing striations form between the point of 
application of the RF discharge and the anode. These are similar to 
those found at the cathode end of the positive column. Both the high 
voltage pulses and pulse modulation of the current also cause disturbance 
of the discharge over large areas. 

There has been no reported attempt to excite waves of stratifica- 
tion by magnetic fields, by the use of radiation, or by acoustic waves. 

In the acoustic work which has been performed in dc discharges, Suber- 
tova has reported creating standing striations by the use of acoustic 
waves, 1^80 j Carretta and Moore observed acoustic waves which they 
attributed to moving striations although these were later shown 

to originate in the cathode region of the discharge, and attributed to 
changes in gas heating as discharge current was changed. Saxton 

investigated the influence of acoustic waves on the positive column and 
found that the arrplitude of the electron density variation was propor- 
tional to the neutral gas density variation. j Wojaczek observed 

that a pulse discharge at one end of the main discharge caused acoustic 
waves which modulated the intensity of the light emitted from the posi- 
tive column, 

With the exception of a brief explanation by Saxton, there has been 
no analysis of the acoustic interactions by the authors of these works. 
Several other experiments involving acoustic waves have been performed, 
but none have involved dc discharges. Ingard has theoretically analyzed 
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the possibilities of producing acoustic waves from interaction with 
charged particles, but has given little consideration to the inverse 
problem. ] 

The way is clearly open for an extensive investigation of the effect 
of acoustic waves on the positive column of a dc discharge. 

In this work, a pulse discharge is used to excite acoustic waves 
of finite extent. Amplitude of these waves is measured by a thermionic 
gas-filled diode, and interaction in the positive column of the dis- 
charge is observed with a photomultiplier tube and displayed on time- 
space photographs. ^14, 76j 

A theoretical explanation of the interaction of acoustic waves in 
the positive column is given. 
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CHAPTER II 



EXPERIMENTAL APPARATUS 

Figure 2-1 shows the arrangement of equipment used in this experi- 
ment. Acoustic waves were excited by a discharge operated between 
electrodes El and E2 , and observed in the dc discharge between electrodes 
E3 and E4. The exciting discharges were produced by high voltage pulses 
in most cases. Square wave modulation of the discharge current was used 
to compare the effect of propagating a compression or a rarefaction 
phase as the leading component of the wave. Hollow cylinders with dia- 
meter 18 mm were employed to take advantage of the higher efficiency and 
higher current densities which characterize the hollow cathode effect. 

^62, 63, 78, 79 j Cross-sectional area is small, giving minimum inter- 
ference with the propagation of acoustic waves. Electrode El was always 
used as the cathode of the exciting discharge; polarity of the electrodes 
in the observing discharge could be reversed easily to study effect of 
direction of propagation in the discharge. 

Power for the observing discharge was provided by a 0-1 kv supply, 
the working load was a pentode with indirectly heated cathode (RCA 807). 
This was chosen to give ripple-free operation of discharge current. 
Variable resistors were used in the cathode circuit, allowing the grid 
to be biased negatively by current flow and the operating current to be 
selected by this simple change of resistance. A pentode was selected as 
the discharge load since this has been shown to increase the range of 
current in which the discharge is striat ion-free, Figure 2-2 shows 

the current and pressure boundary of striation-f ree operation. A curve 
calculated from the empirical relations in also shown. For a 

given pressure, both the upper and lower current boundaries of this region 
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Figure 2-1. Experimental arrangement, showing Z-modulation and diode circuits. 





Figure 2-2. Striation boundaries in a 2.5 cm diameter, 40 cm 
long Neon discharge. 



16 



were found to be larger than the values given by the empirical relations. 
Formulae for these relations are given in Appendix II. 

For observing the interaction in the discharge, the time-space dis- 
play method described by Stirand, et . al., was employed (Figure 2-1). ^76 
In this technique, the oscilloscope horizontal sweep is triggered by the 
exciting discharge, vertical position of the trace corresponds to posi- 
tion of the photomultiplier along the column, and the intensity of the 
beam is modulated by the photomultiplier signal. With the carriage at a 
given position, variation of intensity of the trace indicates a time- 
varying phenomenon. When the carriage is moved to a new position, the 
vertical position of the trace is changed. These variations occur at 
different times for phenomena which propagate with finite velocities. A 
photograph taken with the shutter open while the carriage is moved slowly 
along the column, records the development in time and space of time- 
varying phenomena. Figure 2-3 is a time-space photograph showing the 
propagation of an acoustic wave from anode to cathode in the discharge 
tube. The anode is at the upper limit of the trace, the cathode at the 
lower limit and time increases to the right. Light intensity variations 
occur earlier at the anode end of the discharge because of the anode to 
cathode propagation. Photomultiplier output is at negative potential, 
causing photographs taken in this manner to show light intensity varia- 
tions in an inverted manner, i.e. , increases in light intensity in the 
discharge correspond to darker areas in the photograph. The unmodulated 
level of illumination of the oscilloscope screen was chosen such that 
very small variational signals could be observed. Frequency response of 
the intensity modulation circuit caused distortion of light intensity 
changes longer than about five milliseconds. 
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t ime 

Figure 2 3. iimc- space ph(;tograph showing propagation of 
an acoust ic wave ( rom anode to cathode in a 
2.7 toT'r, 0.5 ma Neon discharge. Time scale 
is 0.5 ms e c /d i V L s io n , 
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Density variations in the acoustic waves were measured by recording 



changes in the conductivity of a gas-filled thermionic diode. This 
technique, first noted by Dayton, et . al. based upon mobility- 

limited operation of the diode according to the equation 

i = KV^. (2) 



Here K involves the mobility which is a function of neutral particle 
density. This device is calibrated by recording current as a function 
of pressure with voltage fixed. In this case, the density changes are 
directly related to pressure changes through the perfect gas law, as 
temperature is also held constant. The functional dependence of current 
upon pressure can therefore be expressed as 



i 




(3) 



where the value of m is obtained from the curves of figure 2-4. Fraction- 
al density variations result in fractional current changes, which may be 
read as voltage variations across a resistor R: 



= -di ^ -dV 
mi^ m.R.io 

From the data in figure 2-4, for a voltage of 20 volts m is 0.286, 
giving 



d>> 

J> 



3.47 







(4) 



(5) 



with R = 1000 ohms, and i^ in ma. 

In the acoustic waves, compressions and rarefactions are adiabatic, 
with pressure and density related through 




const. 



( 6 ) 
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10 volts 




4 fill f « I ^ ' 
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Pressure - torr 



Figure 2-4. Diode current (multiplied by 2) as a 

function of pressure for constant plate 
VO 1 tage . 
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Thus 



• / ( 7 ) 

P 

and the sensitivity of the device at 1.0 torr, with 20 volts from anode 
to cathode is 

S = — — — = 0.24 — — . (8) 

dp .xC(bar 

Dayton *s method gives a sensitivity of 0.79 mV^bar in nitrogen at 1.0 
torr. This can be attributed to the larger total collision cross-section 
of nitrogen. His calculation may be in error, however, as he seems to 
have assumed an isothermal relation between density and pressure instead 
of an adiabatic relation. 

Experiments showed that operation of the neon-filled diode in the 
experimental pressure range was described by the equation 

i = KV^*^ . (9) 

This indicates that the diode is nearly space-charge limited rather 
than mobility limited. The curve of figure 2-4 indicates that the varia- 
tion of conductivity with pressure is sufficient to allow its use as a 
detector of density variation. 

Acoustic wave amplitudes are expressed in fractional variations of 
neutral gas density, as it is these variations which will affect the 
properties of the plasma. All measurements of acoustic wave amplitudes 
were made with the observing discharge extinguished. 

Figure 2-5 shows the diode signals for acoustic waves produced by 
the high voltage pulse discharge and by square wave modulation of the 
exciting discharge current by the pentode. Density increases are indi- 
cated by a negative change in voltage. In this figure it is evident 
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5 msec/division 

Figure 2-5. (a) Acoustic wave produced by 2 kv on plate of 

high voltage pulser. Density variation is 0.27 per cent, 
(b) Acoustic wave produced by square wave modulation of 
discharge current. Peak density variation is 0.094 per 
cent, difference in static levels is 0.04 per cent. 
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that the high-voltage pulses excite waves of much larger amplitude and 
shorter rise time. Maximum density variation is 0.27 per cent for the 
larger amplitude acoustic wave, and 0.09 per cent for the smaller am- 
plitude wave. A static density variation of 0.04 per cent is also 
created by square wave modulation of the discharge current. 

Translation of the gas along the discharge tube axis by a sinusoidal 
acoustic wave is given by 







C*So*Sint4^t dt 



CS 

= (1 - cosuit) , (10) 

where 

C = sonic velocity (456 m/sec @ 24^C) 

S = condensation (dj>/p ) . 

For the high voltage pulser, a frequency of 1 kHz approximates the com- 
pression phase, and the maximum translation is 1.23 mm when the density 
variation is 0.27 per cent. In a later section, this is shown to be 
important. 

Firing of the high voltage pulser causes an electrostatic inter- 
action between the two discharges. This may result from capacitive 
coupling, or from the change in field configurations near electrodes 
as the pulses are initiated. In either case the mechanism is unclear, 
causing an increase in the signal detected in one instance and a decrease 
in another (figure 2-6). 

Figure 2-6 shows the form of the acoustic waves detected by the 
diode for pulser voltages of 2 and 4 kv. This figure shows that the 
higher voltage pulse produces a larger electrostatic interaction, a 
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a 



b 

c 



d 




Figure 2-6. Traces (a) and (d) are 1001 Hz square waves used 
as base lines for calculation of acoustic wave velocity. 

Nominal time base is 0.5 msec/division. Traces (b) and (c) 
are the acoustic waves observed by the diode for pulser voltages 
, of 2 and 4 kilovolts. Density increase is indicated by a down- 
ward deflection of the beam. Deflections prior to 2.5 msec are 
due to electrostatic interaction. The velocity of the acoustic 
wave in (b) is 440 m/sec, and for (c) is 448 m/sec, if acoustic 
waves are assumed to arise at the moment of maximum light in- 
tensity in the exciting discharge. Density variation in (b) 
is 0.34 per cent, and in (c) is 0.71 per cent. 
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larger amplitude acoustic wave, and an acoustic wave which arrives earlier. 
Propagation velocities of these two waves are 440 m/sec and 456 m/sec 
respectively, if the waves are assumed to originate at the time of break- 
down of the gas. The latter velocity corresponds to the sonic velocity 
of small amplitude waves for a neutral gas temperature of 24^C. This 
apparent difference of propagation velocities cannot be resolved without 
a second diode in the tube or some other device to allow time of flight 
to be measured within the tube. If both waves are assumed to propagate 
with the same velocity, 75 microseconds are required for the formation 
of the acoustic wave by the smaller amplitude pulse. 

Damping of acoustic waves by collisions and heat transfer to the 
walls predominates over collision and heat conductivity losses in the 

volume for the pressure range of this experiment. ^25 j The damping 

—2 —1 

coefficient is of the order of 2 x 10 cm . 
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CHAPTER III 



ACOUSTIC WAVE INTERACTIONS 
A. Translation of the Gas 

Acoustic waves were propagated across a short discharge in another 
tube in order to show the translation of the gas caused by these waves. 
In figure 3-1, static light intensity in relative units is plotted as a 
function of distance from the centerline of the discharge. The shape of 
this curve is approximately Gaussian. Arrows indicate photomultiplier 
locations for traces (b) and (c) of figure 3-2. 

Trace (a) of figure 3-2 shows the current in the exciting discharge 
as a function of time. The trace begins with a decrease of current, 
which causes a rarefaction wave to propagate due to the removal of a 
heat source. In trace (b), the photomultiplier is located toward the 
exciting discharge from the centerline of the observing discharge. This 
causes a brighter area to be moved in front of the slit as the acoustic 
wave arrives, while trace (c) shows opposite behavior of the light in- 
tensity as a darker area is moved in front of the photomultiplier. 
Effects are reversed in both traces as current in the exciting discharge 
is increased and a compression wave is propagated. 

When the acoustic waves were propagated along the axis of the dis- 
charge tube, a similar effect was noted in the standing striations of 
cathode regions. However the light intensity was nearly square wave 
modulated when E4 was the cathode and the exciting discharge current was 
modulated by a square wave. The behavior with E3 as cathode was similar 
to that observed for propagation of the acoustic waves transverse to the 
discharge axis, and did not show this square wave modulation. Figures 
3-3 and 3-4 show the waveforms observed when each of these electrodes 
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Figure 3 - 2 , (a) Current variation in the exc'iting dis- 
charge as a liinct. ion of time; 50 ina/d i vm‘.s ion . (b) 

Light intensity variation on the side of the positive 
colunin of the observing discharge toward tiu' exciling 
discharge. (c) Light intensity variation on the side 
of the positive column of the observing discharge^ away 
from the exciting discharge. Current in the observing 
discharge =3.5 ma . Time base = 2 msec/division. 



?S 




50 msec 




20 msec 



2 msec 



Figure 3-3. Light intensity variations on the cathode 
(exciting discharge) side of the first standing stria- 
tion in the positive column of tube II. E3 is cathode, 
0.7 torr, 16 ma. Time per division is shown to the 
right of each trace. Rarefaction wave arrives first 
in each case. All vertical scales are identical. 

Tube II is the tube used for acoustic experiments. 
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a 




Figure 3-4. Light intensity variations in the head of the 
positive column as a function of time. Time scale: 5 
msec/division. (a) Current variation in exciting dis- 
charge, 20 ma/divis ion. (b ) Anode side of first standing 
striation. (c) Cathode side of first standing striation. 

(d) Negative glow. E4 is cathode, 0.7 torr, 16 ma. Rare- 
faction wave propagated first. 



30 




is operated as cathode. Figure 3-3 also shows that approximately 150 
msec is required for light intensity to return to the original value. 
Amplitude of the compression phase of the wave is decreased as frequency 
of the modulating square wave is increased. 

No satisfactory explanation has been found for the nearly square 
wave modulation of light intensity in the cathode regions when E4 is 
used as cathode. This occurred only when current of the exciting dis- 
charge was square wave modulated. 

B. Interaction in a Homogeneous Positive Column 

Time-space photographs are extensively used in this section to pre- 
sent the several effects which were observed when acoustic waves were 
propagated through the discharge. For these photographs, electrode E3 
is always located at the upper limit of the trace, electrode E4 at the 
lowerv-. Travel of the photomultiplier is 40 cm in all cases, with inter- 
electrode separation 39 cm, and approximately equal overlap at each end 
of the carriage travel. 

The translation of gas discussed in the preceding section gives 
rise to large variations of light output for regions of the discharge 
having standing striations. A standing striation is always located at 
the head of the positive column, and other standing striations which may 
be present are damped in the direction of the anode. Damping is in- 
creased as current is reduced, leaving the column nearly homogeneous 
over 3/4 of its length for currents below about 2 ma. 

Interaction in the positive column is first examined by using 
acoustic waves generated from square wave modulation of the exciting 
discharge current, and secondly with the larger amplitude acoustic waves 
from the high voltage pulser- In the first case, the general effect of 
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compression and rarefaction waves upon the column is established, while 
with the larger pulse an additional effect is observed. 

Figure 3-5 shows anode to cathode propagation of acoustic waves 
produced by square wave modulation of the exciting discharge. A rare- 
faction wave is propagated first in all experiments with the square wave 
modulation technique, and it is seen that this rarefaction causes an in- 
crease in light emitted by the discharge (recall that increases of light 
intensity are recorded as dark areas in the photograph). The compression 
wave causes a decrease in the intensity of light emitted by the dis- 
charge. Waves of stratification also originate from the cathode region. 
These will be discussed in another section. 

If the average energy of the electrons remained constant during the 
compression, the increased collision frequency (at the higher pressure) 
would cause an increase of ionization rate and an increase of the light 
emitted by the discharge. Since the light intensity decreases in the 
compression phase, the electron energy must decrease. Increased elastic 
collision or decreased electric field may explain this decrease of elect- 
ron temperature. These will be examined in a later section. 

Figure 3-6 shows anode to cathode propagation of an acoustic wave 
produced by the pulse discharge. In this figure, the path of the acoustic 
wave is shown as a bright line preceded by a thin dark line. The bright 
line indicates a decrease of light intensity in the compression phase, 
and is caused by a decrease in the electron temperature. The dark line 
indicates an increase of light intensity at the leading edge of the wave. 
Two mechanisms can cause this increase: 

(a) compression with constant electron temperature (phase lag 
of temperature change) 
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Figure 3-5. Time and space development of light intensity in a 1.3 
torr, 2.6 ma neon discharge resulting from anode to cathode prop- 
agation of rarefaction and compression waves produced by square 
wave modulation of the exciting discharge. Lighter areas in 
photo correspond to darker areas in the discharge. Time scale 
is 2.0 msec/division. 
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Figure 3-6. Time-space development of light intensity in a 2.2 
torr neon discharge at 3.1 ma. Anode to cathode propagation 
of a compressive acoustic disturbance. Time scale is 2 msec/ 
division. Thin dark line preceding bright line in acoustic 
interaction indicates a slight increase of ionization rate. 
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Figure 3-7. Time-space development of light intensity in a 0.82 

torr 35 ma neon discharge. Anode to cathode sound propagation. 
Time scale 1 msec/division. Striation speeds are 140m/sec, 

66 m/sec, 33 m/sec respectively. The fastest striations appear 
to be excited by travel of the acoustic wave along the positive 
column. 
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(b) increase of electric field and electron temperature. 

This phenomenon was not observed in all photographs with the high-voltage 
pulser, and was never seen when square wave modulation of the discharge 
current was used for acoustic wave production. 

Anode to cathode propagation of the acoustic wave causes changes of 
ionization rate in time and space which are similar to those present in 
moving striations. Figure 3-7 shows striations which are apparently 
produced during propagation of the acoustic waves along the column. By 
viewing this photograph along the apparent path of the acoustic wave, 
with the eye near the plane of the page, a series of lines can be seen. 
The acoustic wave arrives at the anode 1.6 msec after the beginning of 
the Sweep. At this time, light intensity is observed to increase (dark 
area in photograph). This increase can be followed, with some inter- 
ruptions, for most of the distance of the column. A second dark area 
parallel to this is also visible at the anode end of the photograph. 

The position at which the third dark line should appear is somewhat ob- 
scured by the bright areas associated with the maximum density portion 
of the acoustic wave. Succeeding striation patterns become clearer. 

Striation production by the acoustic waves in an inhomogeneous 
positive column is considered in the following section. 

C. Interaction in an Inhomogeneous Positive Column 

As the acoustic wave encounters standing striations in the positive 
column, changes of the ionization rate at these points also produce 
moving striations. Light intensity gradients are small for all but the 
first few standing striations, so that translation of the gas will pro- 
duce negligible variation of the photomult ipler signal in that half of 
the positive column toward the anode. Figure 3-8 shows the standing 



36 



Anode 



40 cm 



Cathode 




Static light intensity 



Figure 3-8. Light intensity as a function of position in a 0.82 
torr, 15 ma neon discharge. Standing striation patterns had 
similar form at 10 ma, smaller amplitudes. Arrows show approx- 
imate phases of standing striations at which moving striations 
of figure 3-9 originate. (Small variation as light decreases 
to zero at anode end is due to a layer of sputtered material 
on the glass . ) 
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striation configuration in a 15 ma discharge at a pressure of 0.82 torr. 
Figure 3-9 shows the production of moving striations by the acoustic 
wave as it encounters the standing striations. These decay and merge 
directly into the wave of stratification vjhich is excited from the stand- 
ing striations nearer the cathode. Wavelength of the standing striations 
is twice that of the moving striations. 

Amplification of the wave of stratification increases as discharge 
current is brought near a critical current (at which self-excited stria- 
tions occur). Waves of stratification can be excited by smaller local- 
ized changes of ionization rate in these cases. Figure 3-6 shows pro- 
duction of a fast and a slow wave of stratification from the region near 

C 

the second standing striation, and a second slow wave from the cathode 
region. In both cases it is difficult to fix the point of origin of the 
wave. A third slow wave is visible at the right hand side of the photo- 
graph. This is due to oscillation of tube voltage and current as the 
first two slow waves arrive in the anode region. Figure 3-7 also shows 
waves of stratification originating from more than one region of the 
discharge. 

At lower currents, amplitude of the standing striations is smaller, 
and waves of stratification originate only in the cathode region. Figure 
3-5 shows that waves of stratification originating in the cathode region 
due to compression and rarefaction acoustic waves have nearly opposite 
courses of development of light intensity. The change of electric field 
in the cathode region causes the wave of stratification to be excited. 
Since the compression and rarefaction waves have opposite effects on both 
of these quantities, the waves of stratification which are produced should 
be opposite in their course. Theory and experiment have shown that the 
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Figure 3-9. Time-space development of light intensity in a 0.82 
torr, 10 ma discharge. Anode to cathode propagation of the 
acoustic wave. Time scale is 1 msec/division. Striation 
production by sound wave is evident. 
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Figure 3-10. Time-space development of light intensity in a 0.7 
torr 16 ma neon discharge. Anode to cathode propagation of 
acoustic disturbances. Time scale is 2 msec/division. 
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wave of stratification is an ”odd” (signed) effect. 43, 48 J 

The effect of direction of propagation is illustrated by comparing 
figures 3-10 and 3-11. In figure 3-10, the wave of stratification pro- 
duced in the cathode region following anode to cathode propagation of 
the acoustic wave is highly damped. For cathode to anode propagation of 
the acoustic wave, striations appear along its path. This behavior was 
also noted in a 4.0 torr discharge, where striations associated with the 
fast wave of stratification were excited along the path of the acoustic 
wave (figure 3-12). Group velocity of the wave of stratification isc 
160 m/sec, significantly smaller than the sonic velocity, and the stria- 
tions excited by the acoustic wave are attenuated before being amplified 
again in the wave of stratification. 

In figure 3-13, the acoustic wave is shown to decrease the ampli- 
tude of self-excited fast striations. These striations are associated 
with direct ionization of the neutral gas molecules, and the decrease of 
electron temperature caused by the acoustic wave reduces the amplitude of 
the striations. 
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Time-space development of light intensity in a 0.7 
ma neon discharge. Cathode to anode propagation of 
disturbance. Time scale is 2 msec/division. 
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Time-space development of light intensity in a 4.0 
0.2 ma neon discharge. Cathode to anode sound propa- 



gation. Time scale: 1 msec/division. 
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Figure 3-13. Light intensity variation with time at a distance of 
1.3 cm in front of the anode. Discharge currents (a) 2.2 ma 
(b) 2.1 ma; pressure 0.77 torr. Cathode to anode sound propa- 
gation. Time scale 2 msec/division, light intensity traces to 
the same relative scale. Light intensity decreases as compressive 
acoustic wave arrives. 
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CHAPTER IV 



ANALYSIS OF ACOUSTIC INTERACTIONS IN THE DISCHARGE 



A. Change of Electron Temperature 

The experiments showed that the light emitted by the homogeneous 
positive column decreases during the compression phase of the acoustic 
wave, and that this must be explained on the basis of a change in temper 
ature. In some cases the light emitted by the column increases for a 
short time preceding this decrease. 

Since the ions are nearly in thermal equilibrium with the neutral 
gas, neutral gas density perturbations will tend to produce the same 
fractional changes in ion density. These can be expressed as 



= 



N 



+o 



N. 



no 



( 11 ) 



where 





n . « N . 
J| Jo 



( 12 ) 



o 

Frequencies of the acoustic waves were <^10 Hz, giving a character 

o 

istic time of 10” sec for the ion density perturbation. Ambipolar 
diffusion times are of the order of 



2 

"2^ — = 4 X 10*“^ sec , (13) 

a 

indicating that ion density changes will be removed through ambipolar 
diffusion to the walls. However, Saxton observed ion density variations 
of the form given by equation (11). Diffusion of the charged particle 
density perturbation evidently does not take place. 

The increase of neutral particle density also decreases the mobil- 
ity of the charged particles, and increases electron energy losses due 
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to an increase of the collision rate. These tend to bring about a de- 
crease in electron temperature. Electric field may also be varied by 
the density change, and thus lead to an additional change in electron 
temperature. In most cases, electron temperature is related to the re- 
duced electric field (E/p), so that it is difficult to separate the 
effect of field and pressure variations on the electron temperature for 
this experiment. 

Using the relation between electron temperature, pressure (density) 
and electric field given in ^17 j, it is possible to examine the effects 
of neutral gas density variations. This equation is 



_vr \ 

® ■ 3 yr • 



\ 

(14) 



where 



T^ is in eV 
e 

^ electronic mean free path 

f_ fraction of electron energy lost per collision. 
This equation can be written 



T 



e 



nTT* 1 E 
3 -Cf /t- 



(15) 



Von Engel graphs the variation of f_ as a function of E/p in j^l7j . 
The curve is for a constant neutral gas temperature and is therefore 
equivalent to f_ as a function of E/N^. For pressures and fields in the 
experimental range this dependence can be expressed as 



f 



“-0 



(_E_)b 

P 



(16) 



Variation of electron temperature with neutral gas density can therefore 
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be determined from 



In general 
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In the succeeding paragraphs, it will be shown that dE/^Njj^O. With 
this knowledge, the temperature change with density depends upon b as 
follows 



b < 2 






0 



b = 2 



dTe 

^Nr 



= 0 



b > 2 






> 0 . 



For the case of neon in the experimental range of reduced field (0.8 - 
5.0 v/cm • torr) , b = 1.9. Thus the variation of electron temperature 
with pressure is small. 

Electric field variation can be analyzed on the basis of the work 
of Pekarek and Masek. their method, the electron density varia- 

tion is assumed to follow the ion density variation quasi-statically , 
with its form governed by the continuity equation: 
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( 21 ) 



(NV)= Q. 



When perturbations of the form 

n^ = N - N 



-o 



n^ N 



-o 



e = E - h I « 



are substituted, this equation becomes 



0 = D_ ^ +yu Eq + u N -4^ 

^^2 ° o -p- 

+ N E . 

o o ^2 



( 22 ) 



(2 3) 



(24) 



Electric field is related to charged particle density through Poisson’s 
equation 
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(25) 



Since 
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then 

“n 

and equation (24) can be written in the form 



Using the two-sided Laplace transform defined in ^6lj by 

oo 

h(t) e~®^dt , 

o<o 

the equation for the electric field is 




(26) 
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where e is the operator image of e. Equation (30) is to be compared 
with Pekarek and Masek's equation (10) which is 



e = J_ (— ) 

s 4TT€o 




(31) 



where 



(32) 



The acoustic perturbation therefore gives a result similar to letting 
^ for which the change of electric field is zero. 

If the change of mobility is not taken into account, the compressive 
acoustic perturbation can be handled through their calculation for 

= 100 This causes an electric field change of -2.5 x lO*"^ v/cm 

for a fractional density variation of 0.5 per cent. 

Comparing the terms in the last bracket of equation (20), and 
assuming = 2 v/cm 

= 3.5 X 10^® cm-3, 

we find 



^ E 



1.3 X 10“^^ 



_M_ =: 5.7 X 10"^^ . 
Nn 



Thus the terms are of the same order of magnitude. It has been assumed 
that the perturbation of ion density does not decay due to ambipolar 
diffusion to the walls. If ambipolar diffusion does prevent this field 
change, there is still a temperature variation resulting from the neu- 
tral gas density perturbation. Such a change must be attributed to 
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variation of collision losses as the collision frequency is changed. 

Magnitude of the temperature change can be obtained from equation 
(20). Recasting it in the form 




(33) 



of 

the effects fractional changes of electric field and neutral particle 
density are readily seen. For the conditions given in the preceding 
paragraph 



^ = -9.3 X 10"^ . (34) 

Te 

The predicted temperature variation is therefore of the order of one 

3 

part in 10 . Sensitivity of the time space display method is quite evi- 
dent from this value. 

The short duration increase of light intensity which sometimes pre- 
ceded the longer decrease in'* the compression phase may be due to the 
effect of the finite length of the acoustic waves. This finite length 
wave can be described by 



^ sin(kZ -cut) U(t - ) (35) 

for a wave which is propagating from anode to cathode with velocity c, 
and a compression phase at the leading edge. In this equation 



J = N - N ^ , |€| « N 

n no ’ |3| no 



(36) 



and the wave is of semi-infinite extent. An exponential damping term 
should be included to give a form similar to that shown in figure 2-5. 
The change of ionization rate can be expressed as 



dz _ dz ^ dNn 
dt dNn " ^ 



(37) 
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and from 17 
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where 

Thus 



X = eV^/kT 
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(39) 



and the solution will have a transient which propagates. This transient 
may explain the slight increase of ionization rate sometimes observed 
at the leading edge of the acoustic wave. 

B. Excitation of Acoustic Waves by the Discharge 

Ingard has predicted that acoustic waves can be generated in the 
neutral gas of a discharge through collisions with the charged particles. 
j^27j Working from the linearized equations for continuity of mass, 
momentum and energy of the neutral gas, with charged particle inter- 
actions entering as source terms, he obtains the equation 



p = V--?. 

^^2 c2 at 

In this equation 

c sonic velocity 

= V7kT^ 

p variation in pressure 

/ ratio of specific heats 

H heat energy source term 

F 



(40) 



( 41 ) 



momentum transfer source term. 



This equation shows that acoustic waves can be excited by a change of 
the heat transfer rate or a spatial gradient in momemtum transfer. 

Ingard attributed acoustic disturbances from moving striations to the 
former, while Briand and Delcroix observed acoustic waves which they 
ascribed to momentum transfer from the ions to the neutral gas. ^ j 

Before carrying out a calculation to compare the magnitudes of 
these two effects in the positive column when moving striations are pre- 
sent, it should be pointed out that variation of the power dissipated in 
the cathode fall due to a change of current will cause a much larger 
acoustic wave than any change in the conditions of the positive column. 
Normal cathode fall for neon with a molybdenum electrode is approximately 
110 volts; “the cathode properties of molybdenum and tantalum are 

similar. Potential difference between the first two standing striations 
in the neon positive column is 18.5 volts, This is larger than the 

potential difference between peaks of moving striations. Assuming all 
energy dissipated in each region is consumed in gas heating, a step- 
wise change in discharge current will cause an acoustic wave to origi- 
nate from the cathode region which will be larger by a factor of 5 than 
that which will be produced from the region between the standing stria- 
tions. However, it is possible to have moving striations present in the 
discharge with no variation in discharge current, or to have current 
variation small enough that the acoustic waves produced in the cathode 
region are smaller than those which arise from the moving striations. 

[lO, ll] 

Consider the case of acoustic waves produced by perturbations of 
an infinite, homogeneous dd discharge column. In order to compare the 
magnitudes of the source terms in the equation derived by Ingard, it is 
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assumed that the electric field is 2.0 volts/cm, electron temperature 
is 2.0 eV, and neutral gas pressure is 1.1 torr. Under these conditions, 
the following parameters are applicable ^ 17 j 

N_ = 2.4 X 10^^ m“^ 

= 1.2 m^/volt sec 
= 80 m2/volt sec 



= 1.2 X 10 



-4 



m 



> = 6.6 X 10-^ m . 

As an approximation to the momentum transfer term, it can be as- 
sumed that all momentum gained by a charged particle from the electric 
field in traveling one mean free path along the direction of the 
electric field is given to a neutral particle at the end of that free 
path. This corresponds to the fact that distribution functions of 
charged particles in electric fields can be divided into a part which is 
symmetric (isotropic) and a part which is asymmetric. The average energy 
of the symmetric part is much larger, and for electrons, corresponds to 
the temperature measurable by an electric probe. Average velocity of 
the asymmetric part corresponds to the mobility or drift velocity. As- 
suming that all particles have the drift velocity of the particular 
species, the momentum transfer is given by 



F = n(m^ . v^ . • Vq ^ • 



(43) 



The effective collision frequency for asymmetric momentum transfer is 



= |vj />, 






(44) 



and as the elastic collision cross-section of neon is nearly independent 
of energy, the mean free path is essentially constant. The velocities 
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used in calculating these collision frequencies are mobility velocities 



= + 



E 



(45) 



since only collisions along the direction of drift motion are to be 
considered for momentum transfer. This gives 

■p = n(. ) E . |e| (46) 

Comparison of the magnitudes of ion and electron contributions to this 
equation shows that 






(47) 



and the electrons can be neglected. Note that this is independent of 
electron and ion temperature. Thus the equation of interest is 



F = 2 






AE 

AX 



(48) 



Large changes of electric field are observed in the presence of moving 
striations, and an assumption of 



A E T / 2 

— — = 1 VO It /cm 

A X 

is well within the range found in other experiments. |^12, 67 j With 
these values, 

V- F = 38.4 joules/m^ (49) 

To compare this with the heat energy term resulting from the 
change of electron temperature in moving striations, it is assumed that 
electron temperature is doubled in 20 microseconds, that each electron 
has the energy of the distribution function average and that all electrons 
collide elastically with neutral atoms. The fraction of energy given up 
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by each species in an elastic collision is denoted f^. Relative con- 
tributions of ions and electrons can be approximated by 



H- 



n« • . Vc- 



kT_ 



= 26 • 



(50) 



H+ n+ • f+ • Vc+ • kT+ 

Thus the ion contribution can be geglected, especially since the tempera 
ture of the electrons doubles while that of the ions changes little. 

Ingard gives the following equation to express the heat contribu- 



tion of the electrons 



[27] 
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Assuming the temperature change to be the dominant factor, this equa- 
tion becomes 
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where 

<d'> = 6 • ( / - !)(^)"^<<r) = 5.2 <<r> 

<<^> = 3.95 X 10"^° 

Ng = 2.4 X 10^® m“3 

c =456 m/sec 

= 300 °K. 

Using the values, and the postulated temperature change, 

= 3. 52 X 10^ joules/m^. (5 3) 

c2 c)t 

This is an order of magnitude larger than the contribution from the 

A P 

momentum term. Discharge conditions under which E, and aflkt 

are larger and A Tg is smaller, will change these comparative magnitudes 
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The experiment of Briand and Delcroix was performed with a plasma 
produced by an discharge. The voltage between two grids bounding 
this plasma was varied at a given frequency, and a microphone located 
normal to the plane of these grids. Acoustic waves were detected at 
the frequency of this modulating voltage. Since the rate of heat gen- 
eration is unchanged by the small modulation signal, these acoustic waves 
must be generated by the term V* ’F, i.e., a net transfer of momentum to 
the neutral gas by the ions. Amplitude of these waves was observed to 
increase as the R-F pow^r increased. Charged particle density increases 
as this power increases, and as V* F is linearly dependent upon this 
density, the experiment and theory are in agreement. 
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APPENDIX I 



DISCHARGE PROCESSES 

This appendix is included for the reader who is unfamiliar with 
glow discharges. The discharge is described generally, and then dis- 
cussed from a microphysical point of view. Particular attention is paid 
to processes important to wave of stratification theory. Cobine 
von Engel [i^j and Loeb ^37j all treat the discharge processes in more 
detail, and the stratification work has recently been reviewed by Oleson 
and Cooper 

A. General Description 

Figure 4-1 shows schematically the major areas of a low-pressure, 
cold cathode dc glow discharge. A discharge of this type in any vessel 
has a cathode fall of definite length determined by the gas pressure and 
type, and by the electrode material. This region is essential to the 
production of electrons for the discharge. The positive column is the 
next important region. A glow discharge can be operated without a posi- 
tive column. Whenever the discharge is longer than the cathode fall, 
the positive column fills the remaining space between the negative glow 
and the anode. It is this part of the discharge which is used for plasma 
research when a weakly ionized medium is suitable. Commercial applica- 
tions of the positive column use its light output in "neon"’ signs and 
fluorescent lights. The anode fall region produces ions for the positive 
column, and is not important unless this fall of potential becomes un- 
stable. In this case, anode spot oscillations develop and disturb the 
properties of the positive column. 

Current in the discharge can be varied over more than three decades 
(0.5 to 500 mA) with constant voltage drop across the discharge tube. 
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ure 4-1. Characteristic light and dark areas of a low pressure 
dc glow discharge. Lower sketch shows electric field as a 
function of position in the positive column, with the coordi- 
nate system chosen such that the positive X direction is in 
the direction of the electric field. 
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This current -voltage behavior characterizes the normal glow regime, and 
finds application in voltage regulator tubes. For discharge currents of 
several hundred milliamperes , heating of a cylindrical cathode becomes 
noticeable. Heating of the tube walls also becomes appreciable at these 
higher currents, particularly in narrow tubes. Sputtering of cathode 
material onto the walls of the discharge tube presents a problem, however, 
particularly as pressure is reduced below two torr. 

The light emitted by each section of the discharge is characteristic 
of that region, and increases with increasing discharge current. Light 
from the negative glow in neon is always orange, as is that from the 
anode fall when this is visible. In the positive column, a definite 
change of color from orange to red is observed as either current or 
pressure is increased. A standing striation always heads the positive 
column, and several other spherically-shaped light regions can be ob- 
served in the column at lower pressures within the experimental range. 
These are highly damped toward the anode. Viewing the striations with 
a prism spectrograph shows that lower energy (red) lines are excited 
closest to the cathode with orange and green lines appearing toward the 
anode. 

Thermionic emission of the electrons from a heated tungsten wire or 
heated oxide-coated cathode eliminates the cathode fall region. Other 
surfaces will produce electrons without heating and allow a low potential 
difference to be used in maintaining the discharge. |^24j Cold cathode 
operation was chosen for these experiments for several reasons: (1) no 
external power supply is needed as in the case of the heated cathodes; 

(2) electrodes may be identically constructed and the polarity reversed 
to check for any directional dependence in the tube; (3) simplicity of 
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construction; (4) insensitivity of this cathode to any environment as 
opposed to the oxide coated cathodes. 

Regardless of choice of cathode, the negative glow is always pre- 
sent and it is this glow that gives the discharge its name. Approximately 
equal numbers of positive and negative charges are present, making this 
the first region with plasma characteristics. 

The plasma of the positive column is a slightly ionized gas, with 

—6 —4 

fraction of ionization 10 to 10 . Collisions with neutral atoms 

play a dominant role in this region. Mobility and ambipolar diffusion 
(described later) govern the motion of charged particles in the column, 
and recombination is assumed to be significant only at the walls. A 
small positive gradient of potential in the direction of anode is essen- 
tial for the production of the charged particles which maintain the 
discharge, and it is this positive gradient which names this portion of 
the discharge. 

B. Cathode Fall 

Electron production at the cathode is accomplished primarily through 
acceleration of ions into the metal surface by the high electric field 
present in this region (figure 4-1). As the ejected electrons are accel- 
erated away from the surface by this same field, they excite neutral 
atoms and form the cathode layer. Further acceleration raises electron 
energy above the ionization level and leads to production of more 
electrons, this time in the Crookes dark space. These electrons repeat 
this chain of events, and the multiplication process grows exponentially 
with distance from the cathode. 

If vie let oc denote the number of ionizations per unit length of 
travel of an electron, and assume that the ionization per length is 
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proportional to concentration, then a single electron emitted at the 
cathode will give 



a 

exp(f o< dx) 

^o 



(54) 



electrons at a distance d from the cathode. Assuming also, that 1/y 
ions must impinge upon the cathode to cause emission of a single electron 
in the steady state each electron emitted must create 1/y ion pairs to 
maintain the discharge. Thus 

d 



y ^exp( r c< dx) - 1 j = 

-^o 



and 



i 



dx = ln(l + — • 



(55) 



(56) 



The integral gives the number of ionizing collisions from the cathode 
to d. This can also be approximated by 

d 

I J cx dx = , ( 57 ) 

where is the number of ion pairs per electron per volt of potential. 



and is the normal cathode fall: 




In (1 + ^ ) « 



(58) 



This equation shows that normal cathode fall increases as ionization 
efficiency decreases, and as the efficiency of production of secondary 
electrons decreases. For the case of neon gas and an aluminum cathode, 
the normal cathode fall is 120 volts, and the length of the cathode fall 
at one torr is 0.64 cm. Tantalum has a slightly lower fall of poten 

tial and shorter dark space, as its characteristics are similar to that 
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of molybdenum. j^3j Complete data on tantalum were not available. 

Cathode potential fall can be reduced and current density increased 
by improving the secondary emission coefficient — the term i in equation 
(58). Although ions have been mentioned as producing the electrons, 
ultraviolet photons can also produce electrons through the photoelectric 
effect. This effect is enhanced by using two large closely-spaced 
parallel plates as cathodes, and is called the hollow cathode effect. 

It may also be obtained by using a cylindrical shell of small diameter, 
and it was in this manner that the hollow cathode effect was used in 
these experiments. Investigations which have been reported show the 
importance of separation distance for flat plates, and the variation of 
this distance with the gas and the pressure. j^78, 79 j Although no data 
was available for the cylindrical shells, the flat plate data gave an 
idea of the minimum and maximum diameters which could be used in the 
pressure range of these experiments. 

Magnetic fields placed such that the magnetic lines are parallel 
to the plane of the electrodes have been used in conjunction with the 
hollow cathode construction. j^62, 63 j These fields have the same effect 
as increasing pressure, provided certain limits are observed. Since an 
equivalent pressure increase can be obtained, the effective lower pressure 
limit for a given geometry can be reduced. 

As electrons leave the high field region of the cathode fall, their 
average energy is reduced. The light emitted in the negative glow shows 
this, as higher energy spectral lines appear nearest the cathode. 

Sputtering of cathode material onto the discharge tube walls hinders 
observation of cathode area light intensity. It has a much more serious 
effect in the case of a sealed tube, for the pressure is decreased as gas 
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molecules are shot into the cathode surface or adsorbed on the surface 
formed by the sputtered material. Variation of sputtering rate is 
functionally dependent upon a negative power of the pressure, as both 
ion velocity and sputtered particle free paths are increased in this 
case, 3 j 

Acton and Swift note that this power seems to be approximately -5, 
and on this basis they derive j^sj 

dp _ -d 
dt “ p5 

P = - 6 

A tube losing 5% of its pressure in two hours will lose 20 % of its pres- 
sure in the next four hours. This effect is known as cleanup and can be 
employed as a means of purifying the gas in a discharge tube. All ele- 
ments other than helium have ionization energies which are lower than 
that of neon, and will therefore be ionized preferentially in any region 
of the discharge. The ions move to the cathode region under the influence 
of the electric field, where cleanup takes place, leaving the body gas 
free of impurities. 

DruyvestB 3 m and Penning note that the Faraday dark space contains a 
**red aureole'* in very pure neon, This is believed to be due to the 

excitation of metastable atoms which have diffused into this region from 
the positive column or the negative glow. It was used as a check on 
purity of the discharge gas in all experiments. A mass spectrometer gives 
the only certain test of the constituent gases of the discharge, and is 
impractical for ordinary discharge work. 

C. Positive Column 

Macroscopic properties of the plasma of the positive column such as 
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the complex conductivity and the spectrum of the emitted light are deter- 
mined by the microscopic processes governing thermal equilibrium, ioniza- 
tion, mobility and diffusion. The development of moving striations is 
dependent upon each of these factors. 

In this section, after a brief description of plasma and the Debye 
length, thermal equilibrium is discussed. Elastic, inelastic and electron- 
electron interactions are included since each of these is important in 
determining the electron energy distribution function. The form of this 
function determines vThether ionization will be predominantly direct or 
stepwise, each case giving rise to moving striations with different 
characteristics . 

The diffusion of charged particles, formation of space charges and 
the effect of electric field changes upon the temperature of electrons 
are all processes which have been studied as a result of theoretical 
papers on moving striations. 

1. Plasma 

1 

Before proceeding with this description of the discharge, it 
should be pointed out that the weakly ionized gas of the positive column 
is a plasma. Irving Langmuir first applied this term to the positive 
columns of discharges with which he was working in order to convey the 
idea that positive and negative charge densities were approximately 
equal. Poisson’s equation 

V • E = 4 TT q^ N+ - N_ j (esu) (59) 

states that large electric fields will be set up whenever N+ / N_, since 
the values of are of the order of 10 cm~'^. The mobility of the 

charged particles will cause this imbalance to be adjusted. The restor- 
ing motion will also cause oscillations, but characteristic frequencies 



64 



for oscillations of both ions and electrons are far above the range of 
this experiment. This internal struggle for neutrality gives rise to a 
characteristic distance for the decay of an artificially introduced 
potential such as at a probe or the tube wall. Known as the Debye length 
or Debye shielding distance, it is given by 




where 

k = Boltzmann’s constant 

= electron temperature, degrees Kelvin 
n^ = electron density, cm*^ 

= electronic charge, esu , 

It must be added that all extraneous fields are not eliminated, 
and the plasma does not succeed in eliminating local electric fields at 

f 

all times. In fact, the question of how the positive column supports 
a macroscopic polarization such as that which accompanies moving stria- 
tions is the very question which has perplexed investigators for years. 
This will be clarified in the section on wave of stratification theory. 

2. Thermal Equilibrium 

The mechanisms involved in attaining thermal equilibrium are 
rather complex, since the electric field acts as an energy source for 
the system, directly affecting both ions and electrons, while these 
charged particles exchange energy among themselves, with each other and 
with the neutral gas. Finally, there is a large energy sink for high- 
energy electrons in the form of inelastic collisions leading to excita- 
tion and ionization. Several attempts have been made to treat this 
problem theoretically, and as each work adds more insight into these 
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fundamental processes, the important aspects of each are given. 

Since the maximum fractional energy which a particle of mass m 2 

can lose in an elastic collision with a particle of mass m^ is 
4mn mn 

^ ^ an electron will lose a small fraction of its energy in 

(m^ + m2) ^ 

colliding elastically with a neutral neon atom - approximately 1 / 37,000 
or 2.7 X 10 “^ - while an ion will lose up to 1/2 its energy in colliding 
with a neutral. This allows the electrons to retain more energy gained 
from the electric field than the ions. Typical average electron ener- 
gies are of the order of two electron volts (eV), equivalent to a kinetic 
temperature of about 23,000 degrees Kelvin. Ions also gain energy from 
the electric field, but ion temperatures differ little from room temper- 
ature due to the efficient transfer of energy to the background gas. 
Thermal equilibrium does not exist between these three species. The 
interaction between them and among particles of a single type all con- 
tribute to establishing a steady state distribution of particles with 
respect to energy. The average energy of this distribution is then, 
of course, the energy associated with the temperature. 

The distribution of electrons with respect to energy changes as 
electric field, elastic collisions, electron-electron interactions, 
and inelastic collisions with neutrals are taken into consideration. 

In the field-free case, the equilibrium distribution is a Maxwell- 
Bolt zmann function, with the temperature of each component equal to 
that of the walls. When an electric field is added with only elastic 
collisions allowed, and the Coulomb interaction between charged particles 
is ignored, the distribution function is known as a Druyvestyn distri- 
bution. 1^,16 j A comparison of these two forms of the distribution 
function shows that the effect of elastic collisions with the neutrals 
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(Druyvestyn) 



X 



o 



= average electron energy 



is to reduce the number of high energy electrons. These high energy 
electrons are essential to excitation and ionization from the ground 
state. Thus the Maxwell -Boltzmann function will cause more excitations 
and ionizations for a given temperature. 

Druyvestyn also derived the distribution function for the case in 
which inelastic collisions are allowed, and Smit’s calculations based 
on this work showed the decrease in high energy electrons which can be 
expected in helium, This loss is reduced as electric field is 

increased. 

Cahn considered elastic collisions and charged particle interactions, 
j^sj His results showed that an increase of electron density caused a 
change from the Druyvestyn toward the Maxwellian form of the distribu- 
tion. A similar result was obtained by Gurevich. j 

Kagan and Lyagushchenko appear to be the first to include both 
elastic and inelastic collisions and electron-electron interactions in 
a calculation of the distribution function. Their result for a 

pressure of five torr and 200 mA is graphed against the Maxwellian 
function and shows that the number of particles at energies above 9 
electron volts is an order of magnitude smaller than the number for a 
Maxwellian distribution corresponding to the same average energy. 
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In each of these cases it is evident that elastic collisions with 



neutral particles cause the deficiency of high energy particles. Re- 
duction of electric field and electron density both serve to make the 
distribution function short of high energy electrons. Thus non-Maxwel- 
lian electron velocity distributions are expected for pressures and 
currents in the experimental range, with the departure from the Maxwel- 
lian fora becoming larger as pressure is increased. 

The distribution function of electrons in neon has been measured 
by Vorob’eva, Kagan and Milenin over the pressure range 0.06 to 1.6 
torr. At 0.06 torr and 200 mA, the distribution was essentially 

Maxwellian. For a pressure of 1.0 torr and 100 mA, the fraction of 
electrons with energy centered at 1 7 eV is half that corresponding to a 
Maxwellian function . This shows the change of the function brought 
about by a pressure change. Data taken at the same pressure, and 25 mA 
discharge current, show that the distribution decreases more rapidly 
with energy than the function obtained at 100 mA. Electron-electron 
interactions tend to make the distribution function Maxwellian. These 
interactions are reduced as current is decreased, and this explains the 
further departure of the electron distribution function from Maxwellian. 

3. Ionization and Recombination 

Ion pairs are produced in the positive column mainly by in- 
elastic collisions of electrons with neutral atoms which are either in 

the ground state or in a meta-stable state. Photoionization plays a 

o 

minor role, as ultraviolet photons with wavelength less than 2400 A 
are required for ionization from metastable states, and vacuum ultra- 
violet photons are required for ionization from the ground state. 

The general equation for the rate of ionization by electrons can 
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be written 






z = 



2k^rl 



m 



[.fl 



(u)*u*f (u)du 
o 



^ J Qs^(u)*u£^(u)dul . 

® Ue. J 



( 63 ) 



The first integral represents direct ionization (from the ground state) 
while the second represents ionization from excited states. Summation 
is over all excited states, although only the four states with config- 
uration 2p^3s have populations which are significant. Two of these are 
metastable. Other symbols not previously used are 



u 



u 



1 

Qi(u) 



dimensionless electron energy (eV/kT^) 
eV^/kT^; is the ionization energy of neon 
ionization efficiency for electron energy u; 
approximated by N.Q^ = 0.055 (V - 20.6) 
density of excited atoms at the level denoted 
by s 



e(V. . V 3 )AT^ 



N.Qi 



N 



= y.0.055(V - “ Vg); 9 / is determined experi- 

mentally, corresponds to a different cross section 
for the respective states 
= neutral particle density 
The equation for expresses the fact that ionization cross section of 
the electrons increases with increasing electron energy. This is true 
for the energies which electrons can obtain under the conditions of 
these discharges. 

In view of the deficiency of high energy electrons noted in the 
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previous section, direct ionization would be expected to contribute 
fewer ion pairs than stepwise ionization. A calculation of the rela- 
tive importance of the two methods of ionization shows that stepwise 
ionization provides more ions by three orders of magnitude than direct 
ionization. ^^9 j Comparison with a Maxwell-Boltzmann distribution shows 
that the distribution function derived with electron-electron, elastic, 
and inelastic collisions taken into account produces fewer ions from 
direct ionization by four orders of magnitude. These results lead the 
authors to conclude that ”the role of direct ionization is negligibly 
small compared with that of stepwise ionization*’. [^0 j This is a con- 
clusion which must be ccontrasted with experimental findings on moving 
striat ion behavior. 

The change of light emitted by the discharge can be considered pro- 
portional to the perturbation of electron temperature, Equation 

(63) shows that the change of ionization rate can also be considered 
directly proportional to the change of electron temperature for small 
perturbations. These changes of light emitted and ionization rate with 
electron temperature are important to mathematical formulation of the 
striation processes, and to interpretation of experimental results. 

Recombination of electrons and ions in the volume of the positive 
column is negligible compared to the rate at the tube walls. Consider- 
ing the requirements for conservation of energy, linear momentum, angular 
momentum and spin, it is evident that the surface atoms of the tube wall 
provide much better opportunities for these conditions to be satisfied. 

It is thought that ions and electrons arrive at the walls separately and 
the recombination takes place through motion along the walls. j 
Since this seems to have no bearing upon the processes of the discharge 
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involved in these experiments, it will not be considered further. 

4. Mobility and Diffusion 

When electrons and ions move through a neutral gas as the re- 
sult of an electric field, their average drift velocities are given by 



1 




= E 






The mobility constant for 
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/ electrons is 
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e ^_ / m~ 
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These are results of derivations by Langevin, and the ionic mobility 
is quite satisfactory. Loeb improved the electron mobility 

expression, but the formula given above provides values which are satis- 
factory to an order of magnitude, and has the advantage of being more 
easily understood. Using values of mean free path and average electron 
energy given in chapter IV, the following values of mobility can be 
calculated: 

= 1.24 X 10^ (cm^/volt sec ) 

= 80 X 10^ (cm^/volt sec ) 

Note that the electron mobility is much higher than the ion mobility. 
This is due primarily to the lower electronic mass, as may be seen from 
a comparison of the two mobility formulas. 

Mobility is also important in the diffusion processes of the posi- 
tive column. Since most recombination of ions and electrons produced 
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in the column takes place at the walls of the discharge tube, a radial 
concentration gradient exists in the tube. Electrons and ions tend to 
diffuse toward the walls with their motion governed by concentration 
gradients and electric fields according to the equation 



= - V (D^ Np t E . (66) 

The opposite signs arise from the difference in signs of the charges. 
In the steady state, particle flux of each type must be equal, and the 
plasma condition forces N . If the electron temperature across 

the tube is constant as is essentially the case, D+ is no longer a 
function of position, ^lyj The electric field can then be eliminated 
from the equation giving 



V = . 57N 

+ytl N 



(67) 
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Diffusion of this type is called ambipolar diffusion, and if >> 

(as the case in this work), the coefficient of ambipolar diffusion 
can be simplified 17 to 



D 



a 






(69) 



Mobility and diffusion play an important role in forming local space 
charges in the positive column, which then change the electric field, 
electron temperature and in turn the ionization rate. Pekarek has in- 
vestigated this, beginning with continuity equations for ions and 
electrons and the Poisson equation, By assuming perturbations of 

the form 
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N+ = + n+ \n\<^ 
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E = Eq + e 
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these equations can then be written 
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Since ion concentration changes much more slowly then electron concentra- 
tion, the ion concentration is assumed given. Substituting equation (74) 
into (73) results in the differential equation 

A^n_ d^- 

- 0 - - 5 ^ 

= Q * '•''Oo "o "+ 

Under the assumption that Q = 0 and a Dirac S' -function initial ion per- 
turbation, it is shown that the relaxation time for the electron distri- 
bution to be established is given by 



-Z-= ^ 



D. 



For the conditions of this experiment 
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5.0 X 10~^ cm 



1.8 X 10^ cm^/sec 



2. 8 X 10 sec 



D- = 

r. = 

Characteristic frequencies for moving striations are of the order of 
0.5 - 50 kHz. Electron relaxation time is several orders of magnitude 
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shorter than times associated with changes of striation parameters. The 
electron density thus follows that of the ions rapidly enough to consider 
it quasi-static. With this result, the electron density can be written 
in terms of the ion density as 



where 



and 



n^(x, t) = 
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and the space charge is then given by the equation 



J^(x, t) = ^ " I ^ 

Equation (81) shows that the space charge is different from zero when- 
ever n^ ^ 0. Its magnitude depends upon n^ and the plasma parameters 
through equation ( 78 )* 

There is another consequence of the ion perturbation. The space 
charges which are formed create electric fields which perturb the main 
dc field. Since the electric field determines the electron temperature 
through the equation 

~<3^e _ 3 _ 16 

(3x 2 3XT e->2 ^®2) 

e 

where 

= electron energy in eV 

= electronic mean free path 
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K* = mean energy loss of the electron in an 
elastic collision , 

the change of electric field also changes the electron temperature. 
Pekarek analyzed this change as a function of distance by introducing 
perturbations of electric field and electron temperature through the 
relations 

u = + e |0|«Ue (83) 

e I 1 

E = Eq + e I e| « Eq (71) 

The spatial change of electron temperature becomes 



-4|- = a^e - bj_e 



(84) 



if second order terms are neglected. In arriving at this relation, the 

quantity is equated to the energy gained from the electric field 

o 

in a mean free path. The constants a^ and b^ are 
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The importance of equation (8^ is the shift of electron temperature 
maximum from electric field maximum. Pekarek has shown that this spa* 
tial phase difference is essential to amplification of the wave of 
stratification. ^51 j 

Since a change of electron temperature also changes the light 
emitted by the discharge, variations in light emitted along the dis- 
charge can be attributed to temperature gradients, local changes in 
electric field, and ultimately to regions of local space charge. 
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The formation of space charges and electron temperature gradients 
also changes the diffusion processes. During the derivation of the 
ambipolar diffusion coefficient, it was assumed that the electron dif- 
fusion coefficient was constant as electron temperature is nearly con- 
stant across the tube. A temperature gradient along the column means 
that a perturbation of ion density will decay through electron tempera- 
ture diffusion as well as through ambipolar diffusion. In order to 
compare the magnitudes of these terms, Pekarek and Krejci write 

equation (66) in its full form 



_ = (D_) VN_ + (N_) VD_ + E yW_N. 
and denote the ratio of the terms of interest by 

N.( VD.) 



R = 



(87) 



D ( ) 



noting that the diffusion previously discussed occurs when 1. Ex- 

panding the diffusion coefficient about the equilibrium temperature 
and making use of the relations 



D_ = kT_^_ /q^ 
0 = k(T_ - T_q) 



the expression for R becomes 
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where 



/ = 1 * 

a kT_ T_o 



(91) 



Except for the case of elastic collision cross-section increasing 
linearly with electron energy, (making /= 0), ■/ is of the order of one, 
and may be positive or negative. Thus determination of R depends upon 



76 



the quotient V0/ Equation (77) gives the form of the electron den- 
sity curve, and a similar equation describes the electron temperature as 
a function of position 



/ Gq (x - I ) n+ ( I ) d I 
- oo 



(92) 



where 
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Calculating the gradients of 0, n. for a sinusoidal ion density perturba- 
tion with wavelength much larger than gives the result 

R = — ^ y . (94) 

^o 

An earlier paper showed 

bi = 3q^/2 . 

The value of R is thus of order one. This means that diffusion of 
electrons by temperature gradient is of the same order of magnitude as 
diffusion due to a concentration gradient. 

Investigation of these discharge processes has come about as a re- 
sult of research into moving striations. Some are applied to this work 
to explain experimental results. Others, such as the diffusion of 
electrons due to a temperature gradient, have not yet been applied to 
striation wDrk. 
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APPENDIX II 



WAVES OF STRATIFICATION 



This section gives a brief description of waves of stratification 
and the various theories which have been set forth before reporting work 
on irradiation of the discharge by an external source. Parameters of 
acoustically excited waves of stratification are also presented. 

A. Theory of the Wave of Stratification 

Several attempts to describe mathematically the development of the 
wave of stratification on the basis of the microphysical processes of 
the positive column have been set forth. The earliest of these works 
was published by Pekarek in 1957. In this approach, the observed 

structure of the wave is postulated, the changes of electric field and 
space charge are related to this structure, and the production and decay 
of the space charges are included by relaxation times characteristic of 
the processes involved. Each maximum and minimum of light intensity is 
indexed so that the light area numbered 2k is actually the k^^ light 
area (figure 1-1). Perturbations of steady state quantities are denoted 
by 



0 electron temperature 
n charged particle density 
m metastable atom density 
q space charge density, 

and the constants ^ , are introduced to describe the decay of space 
charge to zero and the removal of the metastable perturbation by ioniza- 
t ion. 



With these quantities describing the column, and assuming the light 
intensity changes to follow those of the electron temperature, the 
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periodic changes in the plasma are described by two chains, depending on 
whether the ionization is predominantly stepwise or direct. Both chains 
begin with a change of electron temperature which is assumed to be ac- 
companied by a change in ionization rate: 



I. 
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In chain I, direct ionization is assumed to produce an increase in the 
number of charged particles of both signs, and space charges form due to 
the higher mobility of the electrons. These space charges decay to zero 
with relaxation time *2^ = l/o<, but cause a change in electron tempera- 
ture of opposite sign in the adjacent region on the anode side. In 
chain II, metastable atoms produced by the increase in electron tempera- 
ture are ionized in a characteristic time 7^ = 1/^ . From this stage, 
the processes proceed as in chain I. 

It is convenient to use a coordinate system in which the striations 
are stationary for writing the equations which describe these processes, 
and transform the solution back to the laboratory frame (in which the 
striations are moving) after the solution is obtained. Using this method, 
the following equations can be written: 

Chain I 

-^5^ = A’6y - 0<n^ (95) 



Chain II 



dm y 
dt 



B’©y - mp 



(96) 



= y m^ - o< 
dt ^ 



( 97 ) 
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For both chains, the electron temperature is assumed to follow any change 
of electric field immediately. The equation expressing this is 



^ y it) = - it) = -k’- n^ ^ (t). (98) 

The second equality states that the space charges are formed whenever 
there is a change in charged particle density. Substituting this equa- 
tion into the equations describing chains I and II gives the following 
differential equations : 

dB . 

I (99) 

d^e 

II * ( o< + /? ) • (100) 

These equations are solved under the initial condition of a Dirac delta 
function increase of the electron temperature, and the solution gives the 
time development of electron temperature in each of the numbered regions. 
The time of occurrence of the maximum amplitude of each region is deter- 
mined, as are the times of occurrence of half maximum intensity. Since 
the time of occurrence of the maximum of light intensity is related to 
the propagation velocity of the wave (u) and the striation velocity (v) , 
both of which are measurable quantities, and are related in turn to the 
relaxation times through the striation spacing or wavelength (^), the 
relaxation times can be checked from the experimentally determined values 
for agreement with the hypotheses as to the governing processes in each 
chain. This determination of relaxation time has been a most useful 
result. It can be shown that 
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VThere the relaxation time is related to the plasma processes through 



for cases in which there are two processes. Also predicted by this 
theory are 

(a) oddness of the wave of stratification (initial conditions of 
opposite sign produce waves of stratification with light and <ferk areas 
reversed) 

(b) a spatial decrease in amplitude of the wave group followed by 
an increase when the amplification is larger than 1.0 

(c) direct dependence of wave of stratification amplitude upon the 
amplitude of the initial disturbance 

(d) growth of width of the packet in time according to t^. 

These predictions were verified in a later experimental work, with the 
exception of the width of the packet of the fast wave. The experi- 

mental results in this case agreed with theoretical predictions in form 
although not in magnitude. 

It should be emphasized that this theory takes the layered structure 
as a basis, attributes the light intensity variations to changes in 
electron temperature, and attributes ensuing changes in electron tempera- 
ture to the space charges produced by the higher mobility of the electrons. 
Basic microphysical processes in the positive column are not used to 
predict the response of the plasma to the perturbation of ionization 
rate. 

A later theoretical work uses a simplified form of Poisson’s equa- 
tion and the equation of continuity for ions with only production of ions 
due to changes in the ionization rate being considered. The Debye 
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length is essentially assumed infinite, as the electrons arising from the 
perturbation of the ionization rate are assumed to diffuse immediately 
upon being created. Perturbation relations of the form 



E = Eq + e 



|e| «E 



(71) 



N = Nq + n 



|n| «N 



(70) 



substituted into these equations give the equations 




(esu) 
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where is the rate of change of ionization rate with electric field 



and is the electronic charge. The solution to this system of equa- 
tions under the initial condition of a unit step in ion density in both 
time and space is a Bessel function with the product (xt) ^ as the argu- 
ment. With coordinate axes chosen such that positive X is in the direc- 
tion of the electric field, and axes designated as in the time-space 
diagram of figure 1-1 (negative X in the vertical direction), periodic 
solutions arise in the form of a family of hyperbolae (figure 4-2). If 
very short and very long wavelength striations are excluded from con- 
sideration by drawing the dotted lines shown approximately 15 degrees 
either side of the quadrant bisector, the resulting lines have the gen- 
eral appearance of the wave of stratification. This diagram then shows 
the propagation in the cathode to anode direction, the motion of stria- 
tions in the cathode direction, and convex dispersion of the wave (dis- 
cussed later). 

This simple set of equations thus contains the necessary processes 
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for the propagation of a periodic structure with the properties of the 
wave of stratification, provided short and long wavelengths are eliminated. 
An improved treatment must include these processes. Pekarek lists the 
three important physical processes which cause the wave of stratification 
as 

(a) dependence of the idnization coefficient on the electric field 
in the plasma 

(b) the production of space charges due to different diffusion 
rates of the charged particles 

(c) creation of additional electric fields due to the creation of 
space charges. 

Since perturbation of the electric field causes a perturbation of the 
electron temperature, the ionization rate is dependent upon electron 
energy. Later formulations have used this dependence instead of the de- 
pendence upon electric field. Moreover, the shift of electron tempera- 
ture maximum from the electric field maximum has been shown to be essen- 
tial to the amplification of the wave of stratification. |^52 j If the wave 
is rapidly attenuated, its characteristics cannot be measured, while 
large amplification causes nonlinear behavior of the wave, invalidating 
the assumptions of the linearized theory. Thus the behavior of the 
electron temperature is important. 

Pekarek and Krejci base an extensive theoretical work on the equa- 
tions of continuity of ions and electrons, the equation for electron 
temperature variation in a spatially varying electric field, and Poisson’s 
equation. 54, 56 These are written as 
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Figure 4-2. Solution to simplified system of equations 

describing electric field variations in the positive 
column following unit step initial conditions. 
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Perturbations of the form given in equations (70), (71), and (83) have 
been applied. The terms and — (n+ - n_) represent production 

of particles due to a change in ionization rate and loss due to ambi- 
polar diffusion. The term z^Nq 0 was present in the simple set of 
equation (eqs. (102), (103)) which produced the time-space periodic struc- 
ture of f’igure This set included no damping mechanism. Ruzicka 

and Pekarek carried out a solution to the set of equations given above 
with the term z^N^9 deleted, Th^ system in this case gave only a 

monotonic decay of an aperiodic disturbance <in ion density, indicating 
that the change of ionization rate with electron temperature is a neces- 
sary condition for the excitation of a wave of stratification. 

} Tu / iThe equations are simplified for solution in the following 



manner: 



( 1 ) 



d z2 



is neglected as the changes of ion density due to 



additional electric fields will be larger than the diffusion due to con- 
cent rat ion grad ient 
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(2) and ZqNq 0 are neglected in the electron balance equation 

as the change of electron density has been shown to adjust to the change 
of ion density several orders of magnitude faster than the average life- 
time of the charged particles, Electron density is thus quasi- 

stat ionary. 

(3) the system is transformed to a coordinate system at rest with 
respect to the ions by Z = x + 

(4) the term is neglected in the sum + yu . 

The equations which result are: 
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The's'e^^.equat ions receive no further simplification before a general 
solution is obtained. The processes included in these equations are 

(a) change of ionization rate with electron temperature 

(b) creation of space charges due to different free diffusion 
rates of ions and electrons 

(c) creation of additional electric fields due to creation of 
space charges. 

These were shown to be sufficient for the formation of a periodic time- 
space structure. Also included are terms for 

(d) the displacement of electron temperature with respect to 
electric field (necessary for amplification of the wave of stratification) 
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(e) damping of the striations by ambipolar diffusion 

(f) finite diffusion coefficient of electrons 

(g) finite mobility of electrons 

(h) ionic mobility 

(i) effect of electron density oh ionization rate. 

The time and space development of the ion density in terms of these 
basic processes is obtained by using the 2-sided Laplace transform 
described in equation which results is; 

A roc 

-^^Cx, t)=j K(x -J ).n^(f, t).d| - -i- n+ . (109) 

^ -'-oo ^ 



In this equation, the term K(x - J ) is a sum of three exponential 
functions. The equation in this form is very difficult to solve. By 
assuming that the Debye length is nearly zero, the authors are able to 
make simplifications and finally arrive at the following equation: 



cK 

dt 



(z, t) 



kT_ 






T? 



+ 



f 




b 



1 



D-. 






n+ 



~ Ho j ® 



ai(Z - |) 



U(-Z +1 ) n+ (| , t) d| . 



(110) 



Their analysis of each term on the right hand side of the equation is 
important and can be summarized as follows: 

(a) the first term is a diffusion term which damps the very 
short wavelength striations 

(b) the second term makes amplification possible (note that 
b^ = 0 means full damping) 

(c) the integral gives an oscillatory character to the solu- 
tion and also damps the long wavelength striations (through the term 
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3]^) , therefore 

(d) there are no redundant terms in this equation. All the 
major physical characteristics are present, and any additional terms 
can only modify the solution, hut not alter its basic form. 

Solution of this equation under a Dirac 5-fvmction disturbance of 
ion density at time t = 0 has the form 



where 



n+(Z, t) = — r - 0 >(k)t) + 0(k)tj dk 

- oc 



k = 2 TT/X 



<-0(k) = 



^ is striation wavelength 

Aiok 



k2 + 
kl + a^ 



AlO = ^0^l<^Eo + ®1 

0f(k) = -k^Dg + z'bj^ 
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( 112 ) 

(113) 

(114) 



Equation (IIL) is a Fourier integral over all wave numbers. The term 
(?(k) is present as an amplification (or damping) term, and thus shows 
how the microphysical processes affect striation wavelength. Amplifica- 
tion will be maximum for 



^i?(k) 
6 k 




0 , 



and this gives 



where 






n = ^ 

^10 



(115) 



(116) 
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Real values of exist for 0<)\<1, and when used in equation (l]2) give 



10^1 _ Pa a? J_ ^ 

k2 . a2 ■ n ' ' 



(117) 



Wojaczek has discussed a method of solving an integral of the type 
shown in equation (111). [ j Functions cU(k) and 0(k) are expanded in 
a Maclaurin series about the resonant wave number retaining only 

first order terms. This is valid only for large distances from the point 
of origin of the wave (large Z) , and assumes that there is large damping 
for all other wavelengths. Using this technique, the integral in equa- 
tion (lU) can be evaluated, giving 
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The authors plot the solution as a function of time for a given 
coordinate, noting that the plot has the typical form of the wave of 
stratification (figure 2-3(a)) . Group and phase velocities obtained 
from the customary relations 
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are of opposite sign. The wave has anode-directed group velocity, (due 
to the negative sign) and is thus of the general nature of the wave of 
stratification. By using striation wavelength and frequency obtained 
experimentally, together with an estimated relaxation length (1/a^) for 
electron temperature, the quantities calculated from equations (118, J19) 
and ( 120 ) are found to agree with the observed values of these to with- 
in an order of magnitude. This theory therefore provides a satisfactory 
description of the wave of stratification for pressures of a few torr 
and currents of order of an ampere. 

Since the phase and group velocities are of the same order of magni- 
tude, and ju|4|v|, this theory is inapplicable to the slow and fast 
waves observed at lower currents. It also predicts a monotonic decay of 
light intensity on the cathode side of the initial disturbance. 

Lee, Bletzinger and Garscadden have solved equation (HO) by using 
a Fourier series technique. ^36 j Their results show one significant 
difference, in that a forward wave, propagating from anode to cathode, 
is predicted whenever ^^21T/ai^ This has recently been observed in ex- 
periments on a 35 torr argon discharge. ^20 j 

Wojaczek has made theoretical and experimental investigations of 
small amplitude moving striations and waves of stratification. ^84-86, 
88-90 J. Using conservation equations for electrons and lions, Poisson’s 
equation and the equation for energy balance of the electrons (including 
thermal conductivity), he assumes small perturbations of the variables 
with plane wave time and space dependence. A dispersion relation is ob- 
tained, and with the assumption of a complex wave number, amplification 
as a function of frequency is derived, j Comparison with experiment 

shows good agreement in form for both the amplification and dispersion. 
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A plot of R • COq against R • p^, where 
R = tube radius 

= resonant frequency 
Pq = pressure 

gives a curve which is in good agreement with experiment for small 
values of R • p^. Pupp’s earlier experiments also gave data which fit 
this curve, 

Dispersion is the variation of phase velocity with wavelength (or 
frequency). A dispersion curve is usually constructed with the fre- 
quency (^) as a function of the wave number (k = 2 T^/A) • Figure 4-3 
shows such a curve. Since the phase velocity and group velocity can be 
calculated from the relations 

V = Uf/V, (119) 

u = d oJ/6 k, (120) 

the dispersion curve is useful in determining the propagation character- 
istics of the wave of stratification. Krejci and Pekarek have analyzed 
the combinations of backward and forward waves and their associated dis- 
persion curves. ^34j In their work it is shown that the curvature of the 
striation paths also determines the curvature of the dispersion curve at 
the point of the resonant wavelength. Figure 44 shows the idealized forms 
of the waves of stratification in (a) neon, and (b) a neon-hydrogen 
(^ ^/l) mixture. The ratio of group and phase velocities in each wave 
are approximately the same*, magnitudes are^difffereht'. . Dispersion is 
different in that the lines are curved away from the origin (convex) in 
(a), and toward the origin (concave) in (b) . The portion of the disper- 
sion curve corresponding to each of these cases is shown in figure 4-4. 
This dispersion curve is predicted by the solution to equation (Ul) . In 
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Figure 4-3. Dispersion curves given by 

(a) simplified model with infinite Debye length 

(b ) finite Debye length theories. 

Diagram after Pekarek. [ 60 ] 
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Figure 4-4, Time-space development of light intensity in 
discharges in (a) neon, (b) neon-hydrogen, following 
an impulsive perturbation of ion density at the point 
X = 0, t = 0. 
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most cases, however, the resonant wavelength is such that the convex 
portion of the dispersion curve is observed. Curve (a) of figure 4~3 
is predicted by the solution to equations (]D3) and (104). The assumption 
in the latter case was an infinite Debye length, while in the curve 
designated (b) , a zero Debye length was assumed. Thus the plasma pro- 
cesses have not been formulated well enough to allow information to be 
derived from the dispersion curve. The authors outline a method for 
obtaining data from the photographs and constructing dispersion curves, 
and this both sets a standard and gives a simple method of accomplishing 
this work. 

Pekarek offers an explanation of the concave dispersion in the neon- 
hydrogen mixture, and of the anode-directed striation motion in pure 
hydrogen. ^59, 60 j The basic differential equation (eq. 110) has a quad- 
ratic ionization term, and in the case of a finite Debye length, this 
term may shift the ionization maximum toward the anode from the point 
of maximum electron temperature. Since the ionization maximum is shifted 
toward the anode, electric field will also be increased in this direction, 
causing the temperature to increase, and the chain of events to recur. 

This is the tentative explanation which he offers for anode-directed 
striation motion in hydrogen. For the case of the mixture, contribution 
from the quadratic ionization term increases as striation amplitude in- 
creases, giving a larger anode-directed component to the striation 
velocity at the larger amplitude phase of the striation. Thus, the case 
of concave dispersion is thought to be associated with non-zero Debye 
length. 

B. Striation Boundaries 

An empirical relation for the upper and lower striation boundaries 
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is given by Rutscher and Wojaczek* 
upper 







+ B, 




Using the formulas given for the 



and lower 



" p - B, 



boundary currents, it is possible to make calculations for the tubes 
used in this experiment. (Subscripts 1 and s come from the German words 
**langsam" and ”schnelle'*, meaning slow and fast, respectively. Slow and 
fast striations usually arise at the respective boundaries.) Quantities 
are defined as follows for neon, with p given in torr and tube radius in 
centimeters : 



Radius 


^1 


^1 


X 


1 


8.0 


.1 


2.43 


2 


6.6 


.07 


1.83 


3 


6.4 


.06 


1.59 



A 



s 




B 



s 




ag = .24 
‘X = .63 
bs = .93 
P = 1.1 



For tube I, R = 1.0, and for tube II, R = 2.5. Tube I was used in 
the irradiation experiments, tube II in the acoustic experiments. If we 
assume the following values for tube II, which can be obtained from this 
table by interpolation, then the striation boundaries can be calculated 



= 6.5 ;) = 1.71 

= 0.65 
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Figure 4-5. Boundary currents for striat ion-f ree 
operation as a function of pressure 
in tube I. Diameter: 1 cm, length 30 cm. 
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for each of the tubes. Results are plotted in figures 2-3 and 4-5 as 
dotted lines, and the experimentally determined values are shown as points 
with a curve drawn through to allow comparison with the results of this 
calculation. It is evident that the striation boundaries found in this 
experiment are much different from those predicted by the empirical rela- 
tionships. In all cases the experimental boundaries occur at larger cur- 
rents, and are as much as an order of magnitude larger in the case of the 
lower boundary. With the exception of the one cm tube at the upper 
boundary, all experimental points fall in such a manner that a straight 
line is the most obvious approximation to the current-pressure dependence 
on the logarithmic plot* 

C. Waves of Stratification Excited by Acoustic Waves 

Table A-1 shows the general increase in both velocity of propagation 
of the waves (Up, u^) and speed of the individual striations (Vp, v^) 
as pressure is decreased. Subscripts p and r are commonly used to desig- 
nate the slow and fast waves. The fast wave recorded at 0.82 torr may 
have been an s wave. Wave propagation velocities are accurate only to 
16% to 20% because amplification is large in many cases, and at times 
two waves have probably merged. These factors make it difficult to 
determine the maximvun of the wave group amplitude, and the packets are 
wide enpugh that the method given by Krejci and Pekarek in [34j cannot 
be successfully applied. In cases where the sound wave was propagated in 
the cathode to anode direction, this acted to increase the apparent group 
velocity, and the acoustic wave also gave rise to the striations. Stria- 
tion velocity measurements are accurate to 5 to 10%, and wavelength 
accuracy is about 10%. 

Four general variations in parameters of the wave group can be obtained 
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Table A-1 



Characteristics of Acoustically Excited 
Waves of Stratification 



Pressure 


i 


u 

P 


'"p 




u 

r 




\ 


(torr) 


(mA) 


(m/s) 


(m/s) 


(cm) 


(m/s) 


(m/s) 


(cm) 


4.0 


0.3 


36 


8.2 


2.4 


160 


35 


4.4 


4.0 


0.2 


X 


4.8 


2.0 


140 


28 


4.2 


3.3 


0.6 


64 


6.9 


2.3 


X 


33 


4.5 


2.7 


0.75 


50 


7.0 


2.4 


220 


39 


4.3 


• 

CM 


0.50 


65 


4.9 


2.1 


235 


28 


3.9 


2.5 


0.60 


49 


5.5 


2.1 


X 


32 


4. 1 


2.2 


3.1 


67 


20/17 


2. 8/2. 7 


X 


X 


X 


2.2 


1.6 


57 


10.5 


2.6 


X 


31 


5.1 


2.2 


1.3 


50 


9.3 


2.5 


X 


X 


X 


2.2 


0.7 


61 


5.6 


2.3 


X 


X 


X 


1.9 


1.7 


61 


11 


2.6 


X 


X 


X 


1.6 


1.7 


82 


8.4 


2.6 


X 


X 


X 


1.4 


8.0 


133 


28 


2.9 


X 


X 


X 


1.4 


5.0 


108 


22 


2.9 


X 


X 


X 


1.4 


2.6 


100 


13 


2.7 


X 


X 


X 


1.4 


1.5 


105 


8.5 


2.5 


■ X 


X 


X 


1.3 


2.6 


84 


11 


2.6 


X 


X 


X 


1.0 


9.6 


140 


15 


3.0 


X 


X 


X 


0.93 


4.8 


X 


14 


2.6 


X 


X 


X 


0.82 


33 


220 


70 


3.0 


5000 


250 


5.4 


0.82 


10.9 


160 


36/24 


3.0/2. 6 


X 


X 


X 
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Table A-1 
(continued) 



0.82 


10.9 


140 


24/36 


2 . 7 / 3.0 


X 


X 


X 


0.77 


15 


140 


24 


2.9 


X 


X 


X 


0.77 


10 


220 


30/21 


3 . 1 / 2.8 


X 


X 


X 


0.77 


8.0 


220 


29/31 


3 . 0 / 2. 7 


X 


X 


X 


0.77 


5.0 


220 


21 


3.1 


X 


X 


X 
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from Table A-1: (1) wave propagation velocity (u) decreases with in- 
creasing pressure, (2) striation velocity (v) increases with increasing 
current, (3) striation wavelength decreases with increasing pressure, 

(4) striation wavelength increases with increasing current. There is no 
theory which gives satisfactory values of u and v as functions of current 
and pressure for conditions other than the zero Debye length limit. Early 
theoretical works of Pekarek related these velocities to plasma para- 
meters through the relaxation time : 

._2. L_ 

2 u + V 

In order to compare these predictions of relaxation time with experi- 
mental observations, acknowledge of plasma parameters over the experi- 
mental range is necessary. The measurements of plasma parameters which 
were carried out do not cover a sufficient range of current and pressure 
to allow this comparison to be made. 

D. Irradiation of the Discharge 
1. Experimental Arrangement 

Tube I was filled with 2.6 torr neon, sealed, and mounted 
parallel to a longer 2.5 cm diameter tube, with a centerline separation 
distance of 2.1 cm. A strip of black paper approximately 0.6 cm wide 
was attached to the photomultiplier housing. This encircled tube I to 
prevent light of the illuminating discharge from entering the slit. A 
point 10 cm from the anode was chosen for observing the effects of the 
irradiation upon the light intensity at different values of discharge 
current. Time-space photographs were taken with 27 cm of photomultiplier 
travel. Current in the illuminating discharge was 500 mA. A black 
cloth at times covered the cathode region of the discharge in order to 
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insure that all effects observed were the result of changes in the posi- 
tive column, and in each case the cloth had very little noticeable effect 
on the changes brought about by irradiation. 

Peak value of the light emitted by the striations in the illuminating 

discharge was found to increase as pressure decreased. Light associated 

with transitions ending on metastable states is the important quantity 

in this experiment, not total light emitted. Eure has found that inten- 
o 

sity of the 8115 A line in argon increases with increasing pressure. 

The transition involved in the emission of this line corresponds to that 

o 

involved in the strong 6402 A line of neon, with the principal quantum 
number being four in argon, three in neon. Hence a similar pressure de- 
pendence would be expected, and it was found that illuminating discharges 
had a larger effect upon the discharge in tube I as the pressure was in- 
creased. However, Druyvestyn and Penning quote data which show that the 
fraction of power input to the discharge which is radiated in metastable- 
ended transitions increases as pressure is increased from 0.1 to 1.0 
torr, and decreases as pressure is further increased, Kagan, 

Lyagushchenko and Khakhaev, working in the pressure range 1-30 torr, 
also found a monotonic decrease in the total number of transitions which 
end on metastable states. j^31 j 

The increase of total electrical signal from the multiplier seen in 
figure 4-6 may be attributed to the change of discharge color from red 
to orange as the pressure is reduced. Since the S-4 cathode of the 1P21 
has a quantiim efficiency which is rising rapidly in the color interval 
red to yellow, the same total number of transitions will cause the 
electrical signal to be larger when the discharge is more orange (lower 
pressure) . 
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2. Irradiation Effects on Steady-State Conditions 

Figure 4-6 shows variation of light intensity and potential 
across the discharge tube as a function of discharge current. (A line 
drawing is given in figure 4-7 to assist in explaining this figure.) 

In part (a) of figure 4-6, the discharge is not illuminated. Low fre- 

quency (slow) striations appear as current is increased above two milli- 
amperes, and high frequency (fast) striations appear below 0.5 milli- 
amperes. Potential across the discharge tube increases as current de- 
creases until the fast striations appear at 0.5 mA. At this point the 

potential increases less rapidly, indicating that less power is dissipa- 
ted in the discharge when the fast striations are present. This is in 
accord with Pekarek’s finding that the striated state of the column is 
a more conservative state. j A slight oscillation of the potential 

across the discharge tube in the presence of ^ow striations is also 
evident. 

In part (b) of figure 4-6, the discharge is illuminated by a 500 
mA, 11.4 torr discharge. The light intensity trace shows that both 
striation boundary currents have been increased; the lower current bound- 
ary from 0.5 to 0.85 mA, and the upper current boundary from 2.05 to 
2.15 mA. Striations appearing at the lower boundary have larger ampli- 
tudes than in the unilluminated case, while striations appearing at the 
upper boundary have smaller amplitudes. The striations above the upper 
boundary appear to have maximum amplitudes which are less well defined 
than those in part (a) of this photograph. Moving striations were pre- 
sent in the illuminating discharge, and these were unstable due to anode 
spot oscillations. These oscillations of the illuminating discharge may 
have caused variation of the amplitude of the slow striations. 
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a 



b 



Figure 4 ^^. Light intensity (increasing to right) and 

voltage drop across tube I at 2.6 torr. Horizontal scale: 
0.5 mA/div. Vertical scale: 20v/div. , with the zero of 
potential offset. Upper - no irradiation. Lower - 500 mA 
in irradiation discharge. 
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Zero current 
light intensity 



Current in discharge 
(Scale; 0.5 mA/div) 



Figure 4-7. Light intensity and potential across the discharge 
tube as a function of current. This diagram is a key to 
figure 4-14. 
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Illumination of the discharge by light from a discharge in the same 
gas reduces the population of metastable states by promoting atoms in the 
metastable states to higher energy states from which it is possible to 
make transitions back to the original state or to another state. Since 
the population of the metastable states is reduced, the stepwise ioniza- 
tion rate is reduced, requiring an increase in the rate of direct ion- 
ization. This must be accomplished through an increase in electron 
temperature, which in turn is a function of the electric field. The 
trace of tube potential in figure 4-6 (b) clearly shows this increase 
of electric potential in the striat ion-free and slow striation regions. 

An increase is also present in the current range of the fast striations 
above a current of 0.5 mA. In this range of current the electric poten- 
tial across the discharge tube not only increases less rapidly with de- 
creasing current, but is found to decrease immediately below the stria- 
tion boundary before again increasing as current is further decreased. 

working with this concept of the effect of the illuminating dis- 
charge, the decrease of amplitude of the slow striations and increase of 
the upper current boundary indicate that the slow striations may be 
associated with metastable atoms. Reduced amplitude of these striations 
is explained on the basis of the lower population of the metastable 
levels. Similarly, the increase of fast striation amplitude accompany- 
ing the increase of direct ionization under illumination indicates that 
these striations are associated with direct ionization. Pekarek has 
also stated both of these possibilities from his relaxation time measure- 
ments and from his irradiation experiments. ^47, 4sj 

Total electrical signal from the photomultiplier is higher when the 
discharge is illuminated. A change of color from red to orange is also 
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observed, and the higher quantum efficiency of the photomultiplier at 
shorter wavelengths (discussed earlier) may account for this increase. 

3. Irradiation Effects on Waves of Stratification 

Figure 4-8(a) shows a slow wave of stratification in tube I 
for a discharge current of 2.0 mA. Examination of the striation paths 
from the right hand side of the page (parallel to the paths) reveals 
convex dispersion in the cathode half of the column, and concave disper- 
sion in the anode half of the column. Part (b) of the figure shows the 
slow wave propagating in a discharge with the same current. In this case 
dispersion remains convex throughout the column. Amplification of the 
wave is also decreased by illumination. 

Discussion of dispersion showed that convex dispersion is predicted 
by zero Debye-length theories, and concave dispersion is presently 
thought to be due to a quadratic ionization term which shifts the ion- 
ization maximum in the direction of the anode. It is therefore attribu- 
ted to a non-zero Debye length. Since the Debye length exhibits the 
2 

dependence ^ T^/N^, an increase in and a decrease in both act 

to increase The change of dispersion observed in the photograph of 

figure 4-8 would indicate a decrease in Debye length if this explanation 
were correct. 

Assuming that the irradiation effects can be accounted for on the 
basis of changes in electron density and temperature alone, the follow- 
ing line of reasoning applies: 

a) E is increased 

b) Te is increased 

c) Since j = <^ E , where 

j = current density 
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Figure Waves of stratification in tube I at 2.6 torr 

and 2.0 ma. Time: 0.4 msec/div. , deflection: 27 cm. 

Upper - no Irradiation, lower - 500 mA in illumination 
discharge. 
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cT" = plasma conductivity 
and <=-< Nq • 

A constant j requires a decrease in N^, therefore 
d) T^/Nq is increased, and ^ is increased. 

This, however, is contrary to the expected result. 

Before drawing conclusions as to the validity of theoretical ideas 
relating to dispersion, it must be pointed out that insufficient data 
is available. The irradiation effects must be examined for several cur- 
rents in the striat ion-free region to insure that the phenomenon recorded 
is a general occurrence. Data collected during this work indicate that 
concave dispersion is observed less often than convex dispersion. The 
problem requires further study before a statement can be made. 
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APPENDIX III 



OTHER WORKS ON PLASMA-ACOUSTIC WAVE INTERACTIONS 

Suits evidently published the first paper related to this subject 
in 1935. [sij Using an arc discharge as his main discharge, he excited 
sound waves by a high-voltage pulsed spark discharge at one end of the 
arc column, and detected the acoustic waves by observing the variation 
in electric current drawn in a second spark discharge placed near the 
other end of the column. From the time required for propagation, he 
deduced the propagation velocity and hence the neutral gas temperature 
in the arc. Blum, Dyer and Ingard used this effect of gas pressure on 
current in their construction of a corona discharge for use as an acoustic 
transducer, 6 j The corona is a special case of a glow discharge, and 
current drawn for a constant voltage is a function of mobility, which 
in turn is a function of density. They indicated that this was a 
sensitive microphone, but as nothing has appeared in the literature 
other than their abstracts of papers for meetings of the Physical and 
Acoustical Societies, no quantitative information is available* 

WDjaczek inaugurated the present-day form of this work by reporting 
modulation of the light intensity of the positive column as a result of 
a high-voltage, microsecond pulse, This was discovered during 

investigations of the positive column with regard to waves of stratifi- 
cation. He found that the effect was larger when propagation was in the 
anode to cathode direction. An estimate of the damping showed it to be 
an order of magnitude larger than that predicted by acoustic wave theory. 
Hayess performed this experiment for the express purpose of examining 
the damping coefficient as a function of current and pressure. 25 I 
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Working in xenon, he found that varying the current had little effect 
on the damping coefficient, and an increase in pressure resulted in a 
decrease in the damping coefficient. Wall damping was calculated to be 
the most important influence in the narrow (1.55 cm diameter) tubes used 
in the experiment. This was an order of magnitude larger than damping 
due to . coll is ions and thermal conductivity. Wojaczek had considered 
only the latter two mechanisms in his calculations, and this explains 
the discrepancy noted above. 

Berlande, Goldan and Goldstein found modulation of the light emit- 
ted by an afterglow discharge in a cryogenic plasma. ^4j They attribute 
the acoustic wave production to decay of electron temperature, but 
present no analysis to support this contention. 

Goldstein, Roux, and Dayton introduced a sound wave into a dis- 
charge tube by microwave breakdown and heating of the gas at one end. 
^22j This sound wave was detected by monitoring the variations of 
quantities in three spearate devices: (1) light intensity in an after- 

glow plasma located at some distance from this source, (2) microwave 
power transmitted in a waveguide in which the same afterglow plasma was 
used as a dielectric post, (3) current in a gas-filled thermionic diode. 
The latter is explained in Chapter II. 

Gentle and Ingard used a high-voltage pulsed breakdown for sound 
wave production, and biased Langmuir probes into the ion collection 
region to measure the sonic velocity. In this manner they determined 
gas temperature in a dc discharge maintained in a flowing gas. ^2lj 
Nygaard also used a high-voltage pulse breakdown for sound wave produc- 
tion, observing the arrival of the sound wave by means of a time dis- 
placement of the density-sensitive Tonks-Dattner resonances in an after- 
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glow plasma. 

Ingard and Gentle have shown that density variations in the neutral 
gas will cause corresponding variations in the ion and electron densi- 
ties. Saxton *s work confirmed this. ^72j Ingard ’s later work 

also shows that acoustic waves may be amplified in the discharge. ^26j 

Sodha and Palumbo have theoretically analyzed interaction of 
acoustic and electromagnetic waves in a plasma. They note that 

electromagnetic waves in the plasma can be modulated by acoustic waves 
in such cases. Modulation can also occur upon reflection from a plasma- 
free-space boundary which is affected by an acoustic signal. This can 
become important in upper atmosphere detection of acoustic waves gen- 
erated by explosions and rockets. 

Strickler and Stewart were in the process of investigating the 
effect of square-wave modulation of the discharge in connection with 
moving striations, when they discovered that they had caused the column 
to assume a permanent kink, This was attributed to the excitation 

of a normal longitudinal acoustic mode of the column, and is the most 
striking effect observed so far. 
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